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Abstract 

We examine how initial coherences in open chiral systems affect distinguishabihty of pure versus 
mixed states and purity decay. Interaction between a system and an environment is modeled by 
a continuous position measurement and a two-level approximation is taken for the system. The 
resultant analytical solution is explored for various parameters, with emphasis on the interplay of 
initial coherences of the system and dephasing rate in determining the purity decay and differences 
in the time evolution of pure vs. mixed initial states, process. Implications of the results on several 
fundamental problems are noted. 
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I. INTRODUCTION 



The quantum principle of parity conservation prohibits the presence of the left- and right- 
handed states of chiral molecules as true stationary states, since the Hamiltonian and the 
parity operator commute and thus the eigenfunctions of the parity-invariant Hamiltonian 
should show parity symmetry in the absence of energy degeneracy [1]. However, some 
chiral molecules are observed to be stable for an enormously long time despite their lack 
of parity symmetry. Theoretical efforts to resolve this conundrum date back to Hund's 
proposal in the early days of quantum- mechanics [2]. He formulated a description of the 
stability of chiral molecules in terms of quantum tunneling in a symmetrical double-well 
potential, where exceedingly long tunneling times through the inter-well potential barrier 
were responsible for the stability. However, one can show that the tunneling time estimates 
in observed stable chiral molecules have a wide range, from sub-picosecond to the age of 
the universe. On the other hand, there have been attempts to employ the interaction of 
molecules with the environment, which leads to decoherence, and thus to a suppression of 
coherent tunneling oscillations, to explain the stability of chiral molecules. Intermolecular 
collisions [3] , interaction with photons [4] , and with phonons [5] have been introduced as the 
physical origin of the dephasing process. 

Another aspect of stabilization of chiral molecules is the experimental realization of the 
superpositions of stable chiral states (so-called Schrodinger's cat state). Several proposals 
have been suggested for the preparation and detection of the superposed chiral states utiliz- 
ing phase-controlled ultrashort laser pulses [6, 7], and growing attention has been directed 
towards using lasers to manipulate the molecular chirality [8, 9]. Approaches such as these 
which rely on the interaction between an external electric field and a molecule have become 
more promising with recent advances in laser technology. Considering that the decoherence 
effect is crucial to experimental realizations of coherent superposition of chiral states, it 
is thus of great interest to understand the coherent tunneling dynamics of chiral systems 
coupled to external environments more thoroughly. 

In this Letter, we study the role of initial coherences, in an open chiral system, on the 
distinguishability of initially pure versus initially mixed states as measured by population 
difference between the chiral states, and the vulnerability to decoherence as measured by 
the purity decay. The influence of an environment on the system is modeled by continuous 



position measurements, an approach often used in the studies of the quantum Zeno effect 
[10] and in the studies of dccohcrcncc effects in chemical reactions [11]. In particular, 
Milburn and his coworkers [12, 13] studied a model of continuous position measurements on 
a quantum system in a double-well potential, and within a two-level approximation their 
resultant master equation is identical to ours. However, their goal was the demonstration of 
the quantum Zeno effect in certain parameter regimes and thus an analytical solution was 
obtained only for an initially localized state. In contrast, our goal is the identification of 
the roles of the initial coherences and the subsequent decoherence in the distinguishability 
of pure vs. mixed chiral states and the stabilization (or loss) of the chirality. Accordingly 
we have obtained an analytical solution for an arbitrary initial condition and explored this 
solution for a variety of initial conditions and dephasing rates. Our results provide some 
answers to three fundamental questions arising in open chiral systems: 1) Can we distinguish 
a pure state from a mixture with same populations as those of the pure state? 2) Why are 
some chiral systems stable for a very long time? 3) Why are superpositions of stable chiral 
states not observed experimentally? Distinguishability between the initially pure state and 
the initially mixed state is measured by a population difference between the two chiral states. 
The stability of the system is measured by the purity, which is equal to 1 only for a pure 
state and decays as coherences are lost during an interaction between the system and the 
environment. 

II. CHIRAL SYSTEMS IN THE PRESENCE OF DEPHASING 
A. Formalism (a two-level approximation) 

Consider chiral molecules in a symmetric double-well potential interacting with a bath 
and assume that the effects of collisions between the system and the environment can be 
modeled by the continuous measurement of the system coordinate, say x [14, 15]. Then the 
system density operator p, which is obtained by tracing the total density operator over the 
environment, obeys the master equation [11], 




dp 



I 



(1) 



Here Ho is the system Hamiltonian and F is a coupling strength between the system and the 
bath. The second term in this equation, a multiphcation of F and the double commutator, 
represents the environmental effect on the system, and is expected to destroy coherences 
between the eigenstates of the position measurement operator, x. This loss of coherence, 
i.e. , decoherence, has been argued as the essential ingredient to bridging the gap between 
classical and quantum descriptions of a microscopic system coupled to its environment [16]. 

To invoke a two-level approximation, we assume that for an isolated system the potential 
barrier is very high or equivalently the energy of the system is sufficiently low that the higher 
energy levels are not involved in the dynamics. In addition, since this is an open system, we 
should point out that the position measurement of the system is known to cause its energy 
and its energy width to increase over time [11, 17]. The system energy and energy width 
will eventually become greater than the potential barrier, leading to the breakdown in the 
validity of the two-level approximation. The stronger F, the shorter the time during which 
the two-level approximation will be valid. 

Within the two- level approximation, we need only the two lowest energy eigenstates, |-|-) 
(even parity; ground state) and |— ) (odd parity) with energies and E^. Alternatively, 
we can consider the two locahzed states \L) and \R), found in the left- and right-hand wells, 
respectively. By convention, these states are given by |L) = :^(|+) + |— )) and \R) ~ 
^(1+) — |— )). The free Hamiltonian readily takes the form = -^+|+)(+| + = 
E^{\R){R\ + \L){L\) - d{\L){R\ + \R){L\) where E^ = ^^^^ and 6 = ^^-^ = hu 
with a tunneling frequency ou. For simplicity, choosing the energy origin to correspond 
to Ejn — 0, i.e. , E+ — —E- and using the usual Pauli matrices cr^, cr^, in the \L) 
\R) basis, Hq becomes Hq — —5(Tx- Using the two-level approximation and the evenness 
and oddness of |-|-) and |— ), respectively, we derive the position measurement operator as 
^ = + H.c = XavgO^z ~ ^niin^^zi where H.c denotes the Hermitian conjugate 

and ixmin are the positions of the double-well minima, which for a deep well are very close 
to the average positions of the right and left wave functions (denoted by ixavg)- Finally the 
master equation for the two-level system becomes [12] 

dp I 

= ^[5(7x, P] - l[^z, [^z, P]\, (2) 

with 7 = Fx^ijj. This equation gives the evolution of the tunneling two-level system sub- 



jected to a continual measurement of the system coordinate x as represented by the Pauh 
matrix a^. The second term in this equation represents a measurement of a population 
difference between two states \L) and \R), and leads to decay of the coherences between the 
\L) and \R) states, i.e. , of the off-diagonal matrix elements, plr- Note that here 7 is given 
by "a multiplication of F by x^j^" , which means that the magnitude of 7 is determined by 
the interaction strength (between the system and the bath) and the system potential as well. 
Thus within the two-level approximation, 7 can be increased by increasing the interaction 
strength or by increasing the distance between two wells of the potential. It is intuitive that 
for a system composed of two wavepackets that are initially spatially separated and sub- 
jected to an environment destroying coherences between these two wavepackets, the larger 
the distance between two wavepackets and the larger the interaction strength between the 
system and the environment, the faster coherences between two wavepackets are lost [10]. 

It is also noteworthy that this master equation is identical to that resulting from a two- 
level atom driven by a field undergoing a pure dephasing [18] in quantum optics, as previously 
mentioned in Rcf. [12, 13]. This equivalency emerges where the following replacements are 
made: \L,R) the two energy eigenstates, ou — >■ Rabi-frequency, and 7 — >■ pure dephasing 
rate. 

The density-matrix elements of the states \L) and \R), pll, Prr, and plr obey the 
following set of equations: 

^ = -lmi2co prl), (3) 
-lrai2upR,), (4) 



dt 
dpRL 
dt 



= -4:7 Prl + i<^{PLL - Prr) (5) 



It is convenient to switch to the Bloch- vector representation by defining the real quantities, 
X, Y, and Z as X ^ 2Re (plr), Y = 21m (plr), and Z = prr - pll [19]. (In regard to 
relationship of X, Y, and Z to coherence and localization of the system, it is noteworthy 
that for a pure state, since X^ + Y'^ + Z"^ = 1, the localized state, say \R), corresponds to 
Z"^ = 1, i.e. , X^ -|- = 0, and the strongly delocalized state corresponds to pa 0, i.e. , 



+ 1.) Rewriting Eqs. (3) to (5) in terms of X, Y, and Z: 

f - 

dY 

_ ^ _4^y + 2a;Z, (7) 
^ = -2a;r. (8) 

The resultant equations are of standard form [20]. Thus the study of chiral stabihty is 
reduced to the solution of the 3x3 dynamical system which depend upon two parameters 
u! and 7- Note that the longtime steady-state solution of Eqs.(6)-(8) isX — Y — Z — 
as t — )■ oo. It can be obtained by setting dX/dt = dY/dt = dZ/dt = 0. This stationary 
state, which is reached regardless of initial conditions, system parameters, or dephasing rate, 
represents a fully mixed state, with equal populations in the right and left wells and with 
no remaining coherence between them. 

With an arbitrary initial state [X — Xq,Y — Yq,Z — Zq], the general solution is, for 
a; > 7 (tunneling dominant region), 

X = Xoe-^^'\ (9) 

Y = — —[ujYo{-2cos{st) + — sin(st)} - Zo(s + — ) sin(st)], (10) 
^UJ s s 

Z = e-'^*[Zo{cos(st) + ^ sin(st)} - sin(st)], (11) 

s s 

with s = 2y/u'^ — 7^. The population of the system in either well (prr = ^^,Pll = ^-^) 
oscillates at the modified frequency s, that is determined by relative size of the dephasing 
rate 7 and the so-called tunneling frequency uj, under an exponentially decaying envelope. 
The introduction of decoherence increases the period of the oscillation (since the frequency s 
decreases as 7 increases from zero towards ou) while simultaneously suppressing the amplitude 
via the decay factor e"'"'*/^ in the tunneling dominant region. 

On the other hand, for 7 > a; (dephasing dominant region), the general solution is 
obtained by replacing sin(st) by sinh(st) and cos(st) by cosh(st) with s = 2-^/7^ — cu'^ in 
Eqs. (9)-(ll): 



(12) 



e-27* 4^ 4^2 

Y = — - — [ujYo{-2cosh{st) + -rsmh{st)} - Zo{s + ^)smh{st)], (13) 
^UJ s s 

Z = e-='^*[Zo{cosh(5t) + ^ smh(5t)} - sinh(si)]. (14) 

In this regime for an initially localized state, say \R), the system approaches the long- 
time steady-state monotonically without any oscillations present. Ref.[12] demonstrated the 
quantum Zeno effect in this regime: when the measurement strength exceeds some threshold, 
the system initially prepared in one well remains stuck on a time scale of the measurement 
strength. If we attempt to apply this observation to our open chiral system, then we may say 
that for an initially chiral system the time over which the chirality is preserved will increase 
proportionally when 7 increases beyond a certain threshold (here cu). Indeed, this is what 
we observe in the case with an initially localized state for a dephasing dominant region in 
Sec.II C. Further, considering that previous studies focused on manifestation of the quantum 
Zeno effect for initially localized states, it would be interesting to explore whether and how 
this quantum Zeno effect, namely, the suppression of a transition (tunneling), exists in the 
case of an initially delocahzed state; this will be reported in Sec.II C. 



B. Distinguishability 

In this subsection, we examine how the initial coherences affect, over time and under 
the influence of decoherence, the distinguishability between the initially pure state and the 
mixed state. These comparisons require selection of initial conditions for the pure state 
{Xq,Yq,Zq) and the mixed state {X^,Y^,Z^). The most meaningful comparison arises 
from assigning the same initial population of \L) and \R) in the initial {t = 0) pure and 
mixed states, i.e. Zq = Z^, but choosing the degree of initial coherence to be zero for the 
mixed state (maximum randomness scenario), i.e. = Y^ = 0, while varying Y^ for the 
pure state. Note that is determined by {X^y + {Y^y + {Z^^ = 1, and X^ will not 
affect the distinguishabihty as measured by Z since X and Z are decoupled (see Eqs. (6)- 
(8)). As easily can be seen from Eqs. (11) and (14), the difference between the time-evolution 
of Z of the initially pure state and that of the corresponding mixed state, AZ(= Z^ — Z^) 



is given by: 



AZ{t) = ^e-^^*sin(st) for > 7, (15) 

s 

AZ(t) = -^[e^^'-^-^)* - e-(^'+2^)*] for a; < 7. (16) 

Prom the above results one can discuss the distinguishabihty between the initially pure state 
and the corresponding initially random mixture by measuring the population difference Z 
in the presence/absence of the bath effect destroying the right-left coherences. Note that 
the distinguishabihty is not affected by Xq as expected and that it depends strongly on 
the magnitude of Yq. Clearly, if — 0, AZ{t) — for all time and for all 7. That is, 
if Yq — 0, one can not distinguish the initially pure state from the corresponding random- 
mixed state by measuring Z, regardless of the decoherence effect. However, if Y^ ^ the 
distinguishabihty does exist, although the decoherence effect tends to extinguish it over time 
because the decoherence forces both the initially pure state and the corresponding mixture 
to the maximally randomized stationary state (equal |L) and populations without any 
remaining coherences between the two states at long times). It is noteworthy that our 
result reproduces, in the limit of no bath and no time-evolution of the system, the work 
by Harris and his coworkers [21]. They assumed essentially instantaneous measurements, 
occurring on a time scale short compared to that for tunneling and dephasing, and, in the 
framework of the wavcfunction, showed that no parity sensitive experiment can measure the 
difference between them if the states |L) and are taken to be real for the pure state, 
which corresponds to the case of Y^ — 0. 

The behavior of the distinguishabihty also depends on the relative size of the dephasing 
rate and the tunneling frequency. For the tunneling dominant region (cu > 7), AZ(t) shows 
oscillations that arc exponentially damped with time, and whose amplitude at a given time 
tends to increase with increasing u and Y^ . On the other hand, for the dephasing dominant 
region (a; < 7), |AZ(i)| shows a simpler behavior: it increases up to maximum at t = tmax 
from zero ait — 0, and then decreases to zero as time increases. Here t — t^ax — j~ ln[ _ft^^^ ] 
is obtained such that ^^^|^|t=w = 0. 

Typical behavior of AZ(t) for several values of 7 for a given value of cu (here chosen as 1) 
is shown in Fig. 1. For cu > 7 (tunneling dominant region, here 7 < 1), as 7 increases, the 
period of the oscillation increases and the amplitude decays faster, leading AZ to a faster 



approach toward zero (AZ = means no distinguishability) . For example, for 7 = 0.5 
[upper panel], upon the measurement of population after t = 4, one can not tell whether 
the system was initially in the pure state or in the corresponding mixture. On the other 
hand, for < 7 (dephasing dominant region, here 7 > 1), after hitting the minimum point, 
AZ/Yq just increases monotonically towards 0. It is noteworthy that, as 7 increases, the 
distinguishability disappears more slowly, while the maximum of |AZ/1^^| gets smaller and 
is obtained at earlier time, i.e. , tmax gets smaller. For instance, compare the case 7 — 1.5 
and 7 = 10 in the lower panel. For the case 7 = 1.5, IAZ/Yq^I = 0.275 at t = tmax = 0.42, 
and almost zero at t — 8, indicating no distinguishability for t > 8, while, for the case 
7 = 10, \AZ/Y^\ = 0.0494 at t ^ t^^ ^ 0.14, and still nonzero as 0.0226 (roughly half 
of maximum \AZ/Yq\) at t — 8, indicating the long-time surviving distinguishability. This 
observation suggests that considerably strong dephasing may unexpectedly afford means to 
observe the long-time surviving distinguishability. 

C. Decoherence (purity decay) 

In this subsection we examine how initial coherences affect the decoherence process. In 
order to measure coherence decay, i.e. , decoherence, we use the purity, q defined as [22] 
q = Tr(p^), where Tr denotes a trace over the system of interest. Note that the decoherence 
process does not preserve the purity of the state, that is, the purity of the resultant mixed 
state becomes less than 1, while that of the pure state is 1. We first sketch some features 
of the effects of initial coherences and dephasing rate on the purity decay (decoherence 
process), and then present the results of numerical calculations. 

Noting that for a two- level system purity is given by = ^{X'^ + Y'^ + Z'^ + 1) in the 
Bloch- vector representation and using Eqs. (6)-(8), one can obtain 

| = -47(^^ + n (17) 

From Eq.(17) one can expect that at fixed 7, for short times, purity decay gets faster with 
increasing Xq and/or Yq under the restriction Xq +Yq + Zq < 1 (the equality holds only 
for the pure state). However, for the long time region the effects of Xq and Fq on the 
purity decay are different due to different time-evolutions of X and Y (see Eqs. (9)- (11) 



and Eqs.(12)-(14)). While X exponentially decays to zero (see Eq. (6)), F's behavior over 
time is more complicated due to its coupling to Z as seen in Eqs.(7) and (8). It is notable 
that purity decays invariantly with regard to Xq —Xq regardless of time and dephasing 
rates, but not with regard to Yq ^ —Yq. On the other hand, from Eq.(17) at a given initial 
condition purity is expected to decay faster with increasing dephasing rates 7 for short times 
such that X ^ Xq and Y ^Yq. However, this tendency may change in the long time region, 
depending on conditions. Details of the effects of initial coherences, (Xq, Yq) and dephasing 
on the purity decay are explored below. 

All numerical results are for initially pure states. One could also examine purity decay in 
initially mixed states but the initially pure state maximizes the influence of initial coherences 
and is inherently more interesting in terms of applications. We consider sets of initially pure 
states differing in the degree of R and L admixture. That is, (;{t = 0) = 1 and thus 
Xl + Yq + Zl = 1. Here Zq is chosen as 1 — y^Xg + Yq. Also for all the cases examined u 
is chosen as 1, and all the variables are in dimensionless units. 

Figures 2 and 3 show purity as a function of initial coherences (Xq, Fq) with various 
dephasing rates spanning the tunneling dominant region (7 < 1, (a) and (b)) to the de- 
phasing dominant region (7 > 1, (c) and (d)) at times t — 0.1, and 1, respectively. Several 
interesting observations are in order. First, one can easily see that for a given time and 
dephasing rates, purity is symmetric about Xq = at a fixed Iqi but not symmetric about 
Yo = at a fixed Xq. This implies that as |Xo| increases at a fixed Yq = 0, purity decays 
faster for a given time and dephasing rate, but, however, as lYol increases at a fixed Xq = 0, 
decay of purity may be suppressed. Figs. 3 (b) and (c) capture this feature very well. For 
instance, in Fig. 3 (c) one can easily see that, as Yq increases in a negative direction from 
Yq — with a fixed Xq — 0, purity increases up to some point (i.e. , purity increases with 
increasing llol), and then decreases. Second, for very short times such as i = 0.1 [Fig. 2], 
purity tends to decay faster, as 7 increases and as Xq + Yq increases. This agrees well with 
the prediction of Eq. (17), which is <;(t) ^ 1 — 47(Xq + YQ)t for very short times such that 
X ft! Xq and Y Ki Yq. Also note that this proportionality to Xq + Yq gets distorted for 
the larger 7; the circle in the contour of the purity gets distorted in Figs. 2 (c), (e) and 
(d), (f). This is because increasing 7 accelerates changes of X,Y from Xq^Yq. Third, on 
the other hand, for later times (see Fig. 3), and in the region of larger initial coherences 
(Xq + Yq ft 1), purity still tends to decay faster with increasing 7, but for the smaller initial 



coherences (Xq + Yq ^ 0), purity shows the opposite behavior, with increasing 7 beyond 
cj, i.e. , the purity decay starts to be suppressed. That is, if initially the system is strongly 
localized in one of the two wells (smaller initial coherences), then introducing sufficiently 
large dephasing (considerably larger than tunneling frequency) will suppress the decoher- 
ence process (racemization) , but, however, if initially the system is strongly delocalized in 
two wells (larger initial coherences) , then increasing the dephasing rate could accelerate the 
decoherence process and the purity decay at a rate faster than for the initially well-localized 
system for a given condition. This different behavior of the system in respect to increasing 
7, depending on its initial condition, can be understood by recognizing the role of the envi- 
ronment as destroying the left-right coherences, and thus inhibiting tunneling of the system. 
The tunneling (transition) suppression leads to the stabilization of the chirality for the ini- 
tially localized state (or so-called quantum Zeno effect) and the faster racemization for the 
initially delocahzed state. Implications of this interesting result contribute to answers to 
some of long-standing fundamental problems in chiral systems: 1) why are the chiral states 
of some molecules, once produced, stable for a very long time? 2) why are superpositions of 
stable chiral states not observed? Our result suggests that the system, if it is initially a well- 
localized state, will tend to preserve the chirality in the external environment that destroys 
the left-right coherences, provided that the interaction with the environment dominates the 
tunneling interaction, while, on the other hand, the system, if initially a superposition of 
two chiral states, would become very quickly a racemic mixture in the presence of the very 
same external environment. 



III. SUMMARY 



We have obtained an analytical solution for an open chiral system. The interaction of 
the system with an environment was modeled as a continuous position measurement of the 
system and a two-level approximation was made for the system. The interplay of the initial 
coherences and the decoherence on the distinguishability between the pure/mixed states, 
and the racemization/stabilization of the chirality was clarified. If the two chiral states 
are initially taken to be real for the pure state (i.e. , — 0), no distinguishability was 
shown to exist for all time and the dephasing rate, as measured by a population difference 
between the two chiral states, was essentially zero. Increasing the dephasing rate was shown 



to tend to accelerate the distinguishability decay for ^ 0, and the purity decay over 
time. However, considerably strong dcphasing (beyond a certain threshold) was observed 
to enable the long-time survival of distinguishability, and the chirality stabilization for the 
initially localized state. The results should serve as a prototype for understanding the results 
of a chiral system under more complicated and realistic conditions, e.g. , more energy levels, 
more spatial degrees of freedom, and other types of external environment which give rise to 
chiral or achiral interactions. 

This work may be extended to a chiral system interacting with lasers. Recently there has 
been considerable interest on the preparation/control of molecular chirality [6-9]. However, 
the decoherence effect on these preparation/control scenarios is far from well understood, 
especially when the time scales of system dynamics, and/or environmental effects (energy 
loss, dephasing) are comparable to the interaction time scale between the system and the 
laser field. 
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FIG. 1: AZ{t)/Y^ vs. time for various 7 shown inside the box. Upper panel: 7 = (dephasing- 
free case), 7 = 0.05, 0.5 (tunneling dominant region), and 7 = 1.5 (dephasing dominant region). 
Lower panel: all are in dephasing dominant region. All the variables are in dimensionless units. 



(a) 7=0.025 



(b) 7=0.25 




FIG. 2: Purity <^ vs. initial coherences {Xq,Yq) at time t = 0.1 for various dephasing rates 7: (a) 
7 = 0.025, (b) 7 = 0.25, (c) 7 = 1.25, and (d) 7 = 2.5 ((e) and (f) are 2-dimensional contour plots 
for (c) and (d), respectively). All the variables are in dimensionless units. 



(a) 7=0.025 




FIG. 3: As in Fig. 2 but at time t = 1. 



(f)Y=2.5 




(e) Y=1 .25 




